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Summary 

A viscoplastic model is developed which is applicable to 
anisothermal, cyclic, and multiaxial loading conditions. Three 
internal state variables are used in the model: one to account 
for kinematic effects, and the other two to account for isotropic 
effects. One of the isotropic variables is a measure of yield 
strength, while the other is a measure of limit strength. Each 
internal state variable evolves through a process of competition 
between strain hardening and recovery. There is no explicit 
coupling between dynamic and thermal recovery in any 
evolutionary equation, which is a useful simplification in the 
development of the model. The thermodynamic condition of 
intrinsic dissipation constrains the thermal recovery function 
of the model. 

Application of the model is made to copper, and cyclic 
experiments under isothermal, thermomechanical, and 
nonproportional loading conditions are considered. 
Correlations and predictions of the model are representative 
of observed material behavior. 

Introduction 

The design and development of advanced power systems 
require mathematical models capable of accurately predicting 
short-term plastic strain, long-term creep strain, and the 
interactions between them. Multiaxial, cyclic, and 
anisothermal loading states are normal service conditions, not 
exceptional ones. The main reason for doing inelastic analysis 
is not so much to determine the deformation response of a 
structure per se, but rather to assess its useful service life. 
Results from inelastic analyses are input to life assessing 
schemes. These schemes typically require such quantities as 
stress range, mean stress, inelastic strain range, and 
percentages of time-dependent creep and time-independent 
plasticity that make up the inelastic strain range (ref. 1). It 
is therefore of paramount importance that a viscoplastic model 
be capable of predicting these phenomena with reasonable 
accuracy. 

The theoretical development of viscoplasticity has its origin 
with the works of Stowell (ref. 2), Prager (ref. 3), and Perzyna 
(ref. 4), whose models do not contain evolving internal state 
variables. The field experienced rapid advances in the 1970’s 
when internal state variable models began to appear in the 


works of Geary and Onat (ref. 5), Bodner and Partom 
(ref. 6), Hart (ref. 7), Miller (ref. 8), Ponter and Leckie 
(ref. 9), Chaboche (ref. 10), Krieg et al. (ref. 11), and 
Robinson (ref. 12). Theoretical refinements have continued 
to occur throughout the 1980’s in the models of Stouffer and 
Bodner (ref. 13), Valanis (ref. 14), Walker (ref. 15), Schmidt 
and Miller (ref. 16), Chaboche and Rousselier (ref. 17), Estrin 
and Mecking (ref. 18), Krempl et al. (ref. 19), Lowe and 
Miller (ref. 20), and Anand and Brown (ref. 21). Reviews on 
various aspects of viscoplasticity have been written by Perzyna 
(ref. 22), Walker (ref. 15), Chan et al. (ref. 23), Lemaitre 
and Chaboche (ref. 24), and Swearengen and Holbrook 
(ref, 25). Although this list is by no means complete, it is 
representative of the work done in viscoplasticity, and of the 
attention that it has received. 

This report develops an anisothermal viscoplastic model 
applicable to polycrystalline metals. The concept of internal 
state variables is used, and the associated thermodynamic 
constraint of dissipativity is discussed. Three internal state 
variables are used in the model: one to account for kinematic 
hardening, and two to account for isotropic hardening. The 
evolutionary equation for each internal state variable contains 
two terms: one for hardening, and the other for recovery 
(either thermal or dynamic, but not both). This set of 
evolutionary equations is unique to this work. A prime 
objective in the development of this thermoviscoplastic model 
was to make the model as simple as the physics and intended 
applications allow. The fact that thermal and dynamic recovery 
terms are not coupled in any of the evolutionary equations is 
such a simplification. As a result, the material functions for 
the model are constructed readily using established 
phenomenological relationships. The use of both isotropic and 
kinematic variables is necessary to accurately model the 
anisothermal, cyclic, multiaxial, loading states found in power 
system applications. It might be possible to accurately model 
other applications, such as metal forming (ref. 21), by using 
only one internal state variable. 

A wide range of materials are used in devices that are 
designed to operate in high-temperature environments. The 
designer’s choice of material is often dictated by its ability 
to carry load at high temperatures, or to conduct heat. 
Whenever a design calls for a material of high conductivity, 
copper is often the material of choice. For example, a copper- 
based alloy (NARloy-Z) is used as a liner in the regeneratively 
cooled combustion chamber of the space shuttle main engines. 
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This report describes a thermoviscoplastic model for copper 
with this type of application in mind. As an aid to the reader, 
a list of symbols is provided in the appendix. 

Constitutive Assumptions 

We shall consider the anisothermal behavior of viscoplastic 
materials, such as poly crystalline metals. Each element of the 
material is assumed to be isotropic and to carry no stress in 
its initial virgin condition. However, as the material deforms, 
it may lose all or some of its symmetries. Small strains, 
displacements, and rotations make up the deformation of the 
material. Strain e is decomposed into elastic e e (reversible) 
and inelastic, or plastic, e p (irreversible) contributions 

e = e € + e p (1) 

with no net elastic or inelastic strains occurring in the stress- 
free virgin condition. 

The changing internal structure of a material element is 
characterized by its state ( a , a h T). Stress a is the measure 
of elastic change in internal structure, while the internal state 
variables a, (for i = 1, 2, ..., m) are the measures of inelastic 
change in internal structure. This inelastic change affects the 
future response of the material element. Temperature T may 
also influence the elastic and inelastic aspects of internal 
structure. 

Consider those materials whose behavior under 
thermomechanical loading conditions accepts the following 
general constitutive representation 

1 — 2v 

tr{e) = tr{o) + 3a AT (2) 
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m 

4= I ) H v*j (4) 
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where 


i p = e p (a , ocj> T) (5) 

and 

oti = <*/(*> <*j> T ) ( 6 ) 


such that e p and a, both vanish in the stress-free virgin state. 
(A dot placed over a variable denotes its time rate of change.) 
The deviatoric stress S and deviatoric strain e are defined by 

S-.-I tr(a)I (7) 

and 

e = t- ^ tr(t)l (8) 

where / is the unit tensor. The elastic response is governed 
by linear isotropic thermoelasticity, where E , a, and v are the 
material constants. The inelastic response is considered 
incompressible causing the inelastic strain to be deviatoric. 
Onsager relations (ref. 26) are used to relate the internal state 
variables (or fluxes) a, to their conjugate thermodynamic 
forces (or affinities) A { through a symmetric nonnegative 
matrix of moduli H ip which typically vary with temperature. 
(The coefficient of thermal expansion a is not to be confused 
with the internal state variables a,.) 

The inelastic strain and the internal state variables both 
evolve as functions of state. Since stress rate and temperature 
rate are not state variables, equation (6) implies that sudden 
changes in stress and/or temperature do not alter the internal 
structure of the material. Such constant structure experiments 
provide useful information about the stress and temperature 
dependence of inelastic strain rate. 

Theories of plasticity, creep, and viscoplasticity typically 
use state spaces that incorporate the thermodynamic forces A x 
instead of their conjugate fluxes a,. This is accomplished by 
using equation (4) to rewrite equations (5) and (6) as 

e p = i p (a , Aj, T) (9) 

and 

oti = A p T) ( 10 ) 

and therefore, in accordance with equation (4), 

m m or/ 

a= E H w + E (id 

;=i M=i 

establishes the evolutionary form for the thermodynamic 
forces, independent of the thermodynamic fluxes. The last term 
in this equation accounts for the temperature dependence of 
the Onsager moduli, which can be important in practical 
applications with variable temperature histories. 

This constitutive theory must satisfy the thermodynamic 
constraint of intrinsic material dissipativity (refs. 27 and 28) 
defined by 
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a : 'e p > Afit t (12) 

i= 1 

which is a consequence of the laws of thermodynamics. 1 This 
constraint states that the rate of inelastic work done on a 
material element must equal or exceed the rate of increase in 
energy stored by the material due to a changing internal 
structure. (The notation : signifies a double tensorial 
contraction.) 

A special class of materials (called ft-form (refs. 28 and 29), 
or standard (ref. 27), materials) satisfies dissipativity by 
defining a positive scalar-valued function 


stresses account for isotropic hardening effects. To assure an 
initial condition that is isotropic, 0(0) = 0. As 0 evolves, the 
material develops a flow-induced anisotropy. 

A simplified structure for a viscoplastic model is proposed 
in that there is no coupling between separate thermal and 
dynamic recovery terms in any of the evolutionary equations. 
This facilitates the derivation of material functions presented 
in the next section. The proposed thermoviscoplastic model 
is characterized by the flow equation 

i p = QZ(F) (16) 

^2 


Q = Q(<r, Aj; T) (13) 

that vanishes only in the stress-free virgin state. The rate effects 
of the theory (eqs. (9) and (10)) are then assumed to be 
governed by the generalized normality relationships 


and by the evolutionary equations for the thermodynamic 
fluxes 


0 = ^- 
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and 


OLi = 


do_ 

dAi 


(15) 


which necessarily satisfy dissipativity if the set of all possible 
surfaces (0 = constant) is composed of elements that are 
convex, that are nested, and that contain the origin 
(ff-zlj—O). There are many models in the published 
literature that support an Q-form (e.g., refs. 9, 12, and 30). 
It is not necessary, however, that all models represent fi-form 
materials, since fl-form materials are but a subset of the set 
of all thermodynamically admissible viscoplastic models. An 
example of a model that does not admit an fl-form is presented 
in this report. 


Thermoviscoplastic Model 


and 


X = - [F — G]e P 2 = 
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C- D 
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(19) 


or, equivalently, by the evolutionary equations for the 
thermodynamic forces (see eq. (11)) 


B = H, 



H„ dT 


( 20 ) 



e p 2 -QR(G) 


d m d • 

+ -T 

H d dT 


( 21 ) 


The thermoviscoplastic model presented herein uses three anc j 
internal state variables: they are a deviatoric tensor 0 and two 
positive-valued scalars 8 and X. Their conjugate forces are 
referred to as the (deviatoric) back stress B , the drag stress 
D, and the limit stress L, respectively. The back stress accounts 
for kinematic hardening effects, while the drag and limit 

where 


; fyD _ ^ n LdH ( - 

L = — [F - G e 2 + T 

L H t dT 


( 22 ) 


thermal and intrinsic dissipation are not coupled, since equations (9) and 
(10) define the state space as being independent of the thermal gradient (see 
ref. 27). 



(23) 


3 


and 


G = 


C-D 


E = S-B 








(24) 

(25) 

(26) 
(27) 


Here C, H b , H d , and //< are positive- valued material 
constants, 0 > 0 is the thermal diffusivity, Z = e 2 /0 ^ 0 
is the Zener- Hollomon (ref. 31) parameter, R > 0 is the 
thermal recovery function, F > 0 and G > 0 are flow 
functions, D is the effective stress associated with plastic flow, 
and D 2 and e 2 are the magnitudes of effective stress and 
inelastic strain rate. 

The flow equation, equation (16), is a second invariant 
formulation derived from a potential function (eq. (14), see 
ref. 29) in accordance with the thermodynamic development 
of Rice (ref. 32). First and third invariants are not introduced 
into the model because they result in higher order effects that 
are neglected in accordance with von Mises (ref. 33). This 
flow equation is compatible with the kinematic constructs 
proposed by Prager (ref. 34) in his theory of plasticity. 

The isothermal form of the evolutionary equation for back 
stress, equation (20), was proposed by Armstrong and 
Frederick (ref. 35) for plasticity. Its anisothermal counterpart 
was subsequently used by Chaboche (ref. 10) and Walker (ref. 
15) in their viscoplastic models. This evolutionary equation 
introduces an evanescent strain memory effect (evanescent 
along the plastic strain path) that is consistent with IFyushin’s 
hypothesis (delay-trace hypothesis (ref. 36)) that a material 
response depends not on the whole of the previous strain path, 
but only on the most recent part of it. In this evolutionary 
equation, strain-induced dynamic recovery competes against 
strain hardening. This competitive process causes the back 
stress B to evolve exponentially to an asymptotic limit stress 
L associated with kinematic saturation B 2 sat . A graphic 
illustration of this evanescent concept is given in figure 1 . Here 
the limiting state of back stress defines a hypersurface of radius 
L in state space. Whenever this upper bound is reached, a 
perfectly plastic or creeplike response is attained, provided 
that the material has also isotropically saturated. At kinematic 
saturation, the inelastic strain rate e p becomes coaxial with 
both the deviatoric stress S and the back stress B (in agreement 
with the experimental observations of Phillips (ref. 37)), and 
the nested set of flow surfaces |E 2 = constant] becomes 
stationary until unloading occurs. Otherwise, this set of flow 



Figure I.— State space representation of limiting state of back stress. 

surfaces can translate freely within the limiting surface, as 
governed by the flow and evolutionary equations. 

The evolutionary equation for drag stress, equation (21), 
is consistent with the Bailey-Orowan (refs. 38 and 39) 
hypothesis of thermally induced recovery competing against 
strain hardening. It has been used by the author (ref. 40) in 
its isothermal form. Inclusion of the temperature dependence 
of the hardening (Onsager) modulus H d has not appeared in 
previous viscoplastic models. This evolutionary equation has 
the drag stress harden at the rate H d ‘e 2 /G , which equals 
HJ^C - D)e 2 /L because of equation (24). Therefore, similar 
to the back stress, the drag stress hardens exponentially to an 
asymptotic limit C. This evolutionary equation for drag stress 
therefore has the concept of dynamic recovery built into the 
equation without introducing a separate term for it. However, 
unlike the back stress, the effect of strain hardening can never 
be completely cancelled out by the influence of dynamic 
recovery at the isotropic saturation of state; the influence of 
thermal recovery is also required because C > D for all values 
of drag stress. Notice also that this expression for thermal 
recovery can be cast into a strain-induced dynamic recovery 
format, as witnessed in equation (18). 

The evolutionary equation for limit stress, equation (22), 
is unique to this work. The concept of the limit stress evolving 
as a state variable is not a new idea (as it has been used by 
Marquis (ref. 41) and Walker (ref. 15)), but rather, it is the 
form of the equation that is new. In the models of Marquis 
and Walker, the limit stress evolves exponentially by a 
relationship of the form 


L = H ( 



(28) 


to an asymptotic saturation value of the limit stress L sal . 
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Actually, both authors chose to express the limit stress as a 
function of accumulated inelastic strain by integrating equation 
(28), thereby eliminating the need of explicitly defining the 
limit stress as a state variable. In the model proposed herein, 
the term in the bracket of equation (28) is replaced by a 
difference between the two flow functions, (i.e., D(F — G)/L), 
a result derived from thermodynamic and steady-state con- 
siderations. The thermodynamic constraint of dissipativity, 
which motivated the choice of DIL as the coefficient to the 
difference (F — G), is discussed in the next paragraph. At 
steady state, the flow function F is related to the magnitude 
of deviatoric stress 



by definition through the relation 

^U = | (30) 


experiments with a range of 6 to 1 in stress amplitude (obtained 
at different temperatures) were used. The result was the 
saturated value of limit stress L sat varied with applied stress, 
whereas, the saturated value of drag stress Asat was 
independent of applied stress. The proposed evolutionary 
equation for limit stress, equation (22), is consistent with this 
observation (the stress dependence enters through the flow 
function F, as evident in eq. (19)). Plus, it provides a self- 
consistency to the model in that its flow and evolutionary 
equations reduce to the multiaxial formulation 



at steady state, which is the flow equation used in the classical 
theory of creep (ref. 43). 

For this viscoplastic model to be thermodynamically 
admissible, it must satisfy the constraint of dissipativity given 
in equation (12), which becomes 

a:e^> J?:0 + D6 + ZA (35) 


where ss denotes steady state. The flow function F simplifies to 


Combining this inequality with the constitutive equations given 
in equations (16) to (19) results in the constraint 


F L = 


S 2 — B 2 I S 2 — L I 


R>Z 


Bi 

G + 2 — 

LD 


at steady state, which when combined with equation (30) 
enables one to determine another flow function: that is, 



is the magnitude of back stress. The flow functions Fand G 
are defined to be equal at steady state— a result that must be 
true if equations (19) and (22) are to have meaning (since 
/S = <5 = X = Oat steady state, by definition). The derivation 
of equation (31) makes use of the facts that E 2 = ( S 2 - B 2 ) 
and B 2 = L at steady state (neither of which is true, in 
general, under other conditions). Looking at equation (32), 
one can ask: Which state variable ( L , D, or both) accounts 
for the stress dependence S 2 at steady state? To answer this 
question, a nonlinear, least-squares, optimization algorithm 
(ref. 42) was used to determine the saturated values L sat and 
Z) sat from stable hysteresis loops of stress versus strain, where 
kinematic saturation occurred prior to each reversal. Several 


which bounds thermal recovery from below. This constraint 
will seldom, if ever, be realized in any given analysis, but 
nevertheless, it is a physically mandatory constraint. 2 All the 
parameters in this inequality are nonnegative valued by 
definition; therefore, this inequality is satisfied whenever the 
inner bracketed term exceeds 1 in value, which is typically 
the case. Dissipativity provides a real constraint when choosing 
forms for evolutionary equations. In particular, the hardening 
terms are the ones that cause constraint; whereas, the recovery 
terms are of a thermodynamic benefit. For example, if the 
term D(F — G)/L found in the evolutionary equation for limit 
stress, equation (22), was replaced by the difference (F — G), 
then the constraint of dissipativity would require that 



(37) 


2 In the numerical computations done for this report, this thermoviscoplastic 
model satisfied dissipativity at all states except at the very onset of deformation, 
(i.e., at the initial condition). Whenever the inequality of equation (36) was 
not satisfied, the thermal recovery function R was set equal to the right-hand 
side of equation (36), thereby satisfying thermodynamics. 
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instead of equation (36). The additional term FG[(L/D) - 1] 
is sufficient to cause this inequality to not be satisfied in 
practice (at least that was the author's experience). 

In addition to the constraint of dissipativity given in equation 
(36), thermal recovery R is also constrained by the fact that 
the internal state variable 5 associated with the drag stress D 
cannot dip below a minimum value <5 0 associated with the 
virgin, or annealed, state. 

The constitutive equation for the evolution of back stress 
is isotropic, as presented in equation (20), since a scalar H b 
replaces the general, nonnegative, fourth-order tensor H b for 
the hardening modulus. 3 For anisotropic materials, the back 
stress evolves as 

B = (38) 


forms for its material functions 0, /?, and Z must be deter- 
mined. Physically based phenomenological relations are used 
to meet this need. Conditions of steady-state flow provide the 
form of these functions for the class of materials considered 
(i.e., polycrystalline metals and Class II solid solution alloys 
where solute drag effects are negligible (refs. 8,16, and 20)). 

The temperature dependence of the model is largely con- 
tained within the thermal diffusivity parameter, which is taken 
to be (ref. 45) 



and in accordance with Lee and Zaverl (ref. 44), the flow 
equation becomes 


k p = ez(F) 


ML 




(39) 


where 


F = 

D 


(40) 


and M is a nonnegative fourth-order tensor that is a measure 
of anisotropy. The actual characterizations of H b and M are 
beyond the scope of this paper. Nevertheless, as a first 
approximation it may be reasonable to consider 


H b = H b M (41) 

so that the strength of anisotropy is the same for both tensors. 
The fact that the hardening (Onsager) modulus for the back 
stress is, in general, a fourth-order tensor seems to have eluded 
viscoplastic modelers. 


(42) 


where Q is the activation energy for self-diffusion, T m is the 
absolute melting temperature, and £7 has its usual meaning. 
This relationship represents an apparent activation energy, as 
depicted in figure 2 for creep of copper. The cluster of datum 
points that falls below the main trend of the data (in the range 
of about 0.5 to 0.6 homologous temperature) is attributed to 
dislocation enhanced diffusion (ref. 46). The author’s attempt 
to account for this additional mechanism resulted in a poorer 
correlation of the Zener-Hollomon parameter than that which 
is presented in figure 3. 

Since the steady-state flow of metals can be modeled by a 
power law at the lower stresses and an exponential at the higher 
stresses (ref. 47), one can write the steady-state Zener- 
HoIIomon parameter associated with equation (34) as (ref. 40) 



for S 2 < C 


for S 2 > C 


(43) 


Material Functions 

Before the thermoviscoplastic model of the previous section 
can be used to model a material or a class of materials, specific 


3 Unlike classical elasticity which requires two scalars (e.g., E and v) to 
represent the fourth-order elastic modulus in the isotropic case, the hardening 
modulus H b requires only one (i.e. , H b ) because the back stress is deviatoric. 


where A, C, and n are material constants. (The constant A is 
not to be confused with the thermodynamic forces A ( .) A plot 
of this parameter is presented in figure 3 for copper. This 
Zener-Hollomon parameter is defined so that it is continuous 
and differentiable across the interface S 2 = C (i.e., the 
power-law breakdown stress), and is therefore equivalent to 
the hyperbolic function proposed by Garofalo (ref. 48). 
According to Nix and Gibeling (ref. 49), diffusional recovery 
mechanisms control the creep response in the low-stress, high- 
temperature, power-law domain; whereas, dynamic recovery 
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Figure 2.— Apparent activation energies for creep of copper. Data are as 
presented in reference 46. Solid line represents equation (42) for 
constants in table I. 


mechanisms control the creep response in the high-stress, low- 
temperature, exponential (power-law breakdown) domain. In 
more conventional terms, creep is synonymous with a power- 
law response; whereas, limit plasticity is synonymous with 
an exponential response, where the flow stress appears to be 
strain-rate insensitive over a fairly large range in strain rate. 

One can combine equations (30) and (43) to obtain the 
Zener-Hollomon parameter associated with equation (16); that 
is, 


Z(F) 


(AF n 

for 

F < 1 

[A exp [n(F- 1)] 

for 

F> 1 


(44) 



10° 10 1 10 2 10 5 
TENSILE FLOW STRESS, 0 S 2 * MPa 


Figure 3.— Stress dependence of Zener-Hollomon parameter for copper. 
Data are as presented in reference 40. Curve represents fit of equation 
(43) to data for constants in table I. 


Likewise, noting that F = G and R = Z at steady state (a 
consequence of eqs. (18) and (19)), one obtains the thermal 
recovery function associated with equation (18): that is. 


R(G) 


(AG n 

for 

G < 1 

1a exp[n(G - 1)] 

for 

G > 1 


(45) 


Thus in accordance with Nix and Gibeling, recovery is 
controlled by diffusional mechanisms whenever G < 1, and 
by dynamic mechanisms whenever G > 1 . This thermal 
recovery function is thermodynamically constrained by 
equation (36). It is also constrained to be zero valued in the 
virgin, or annealed, state where 6 - 5 a (if thermodynamically 
admissible). Since the flow functions Fand G are not restricted 
to be steady-state parameters, the material functions given in 
equations (44) and (45) for the Zener-Hollomon parameter Z 
and the thermal recovery function R are assumed to be 
applicable to transient and steady-state conditions alike. 

The thermoviscoplastic model has three elastic constants to 
determine: 

E, a, and v 

all of which can vary with temperature, and it has ten inelastic 
constants to determine: 

A f C, d 0 , H b> H d , H e , \ ot n, Q y and T m 

of which only the three hardening (Onsager) moduli H b , H d , 
and H t can vary with temperature. Two material constants are 
required to characterize thermal diffusivity — Q and T m , three 
to characterize steady state— A, C, and n y three to characterize 
transient behavior—//^, H di and and two to characterize 
the initial condition— b 0 and A 0 . The values of these elastic 
and inelastic constants are given in table I for copper. 
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TABLE I — MATERIAL CONSTANTS FOR COPPER 



a Temperature T is in degrees centigrade; k = 8.314 J/mol K. 


The activation energy for self-diffusion Q and the melting 
temperature T m are handbook data. 

With the thermal diffusivity now characterized, one can 
determine the Zener-Hollomon parameter Z - i p 2 (Q at steady 
state and plot it against its associated flow stress, as in figure 3. 
Each monotonic datum point in this figure represents the 
average of typically four to six creep experiments, usually run 
at different temperatures but with the same flow stress. The 
saturated, fully reversed, cyclic data are in agreement with 
the monotonic creep data. Values for the steady-state material 
constants A , C, and n are obtained by fitting these data to 
equation (43). They are quickly obtained by trial and error. 
The coefficient^ translates the curve along the ordinate. The 
power-law breakdown stress C splits the stress domain into 
power-law and exponential stress-dependent domains. The 
exponent n defines the slope of the curve in the power-law 
domain. 

With the thermal-diffusivity and steady-state material 
constants specified, one can determine the kinematic hardening 
modulus H b from stable, stress-strain hysteresis loops where 
the isotropic state variables have saturated. To accomplish this, 
a nonlinear, least-squares optimization algorithm (ref. 42) was 
used to fit H b and the saturated value D sat to data obtained 
from stable hysteresis loops at 22, 200, 400, and 600 °C. 
Although it is not necessary to use such a technique, it is useful 
and expedient to do so. The experimental and correlated 
hysteresis loops are compared in figure 4. The degree of 
correlation is quite good. 

Four material constants remain to be determined; they are 
the two isotropic hardening moduli H d and H e , and the two 
initial conditions for the isotropic state variables d 0 and \ 0 . 
Once again, a nonlinear, least-squares optimization algorithm 
was used to determine numeric values for material constants. 
This was not necessary, but to determine these constants by 
trial and error would have been a tedious and time-consuming 
chore. Cyclic stress-strain data obtained from the virgin state 
to the cyclically stabilized state at 22, 200, 400, and 600 °C 



STRAIN, € 

(a) Experimental data (ref. 40). 

(b) Theoretical correlations. 

Figure 4 — Saturated, uniaxial, hysteresis loops for copper. Strain rate, 
e, 0.001/s. 
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were used for correlation purposes. An example of 
experimental and correlated stress-strain cycles is shown in 
figure 5. The degree of correlation is quite good. 

The curve depicting the evolution of back stress in figure 5 
provides useful insight into how the three internal state 



(a) Experimental data (ref. 40). 

(b) Theoretical correlation. 

Figure 5.— Uniaxial, cyclic, stress-strain response for copper. 
Temperature, T , 200 °C; strain rate, e, 0.001 Is. 


variables of this model interact with one another. The back 
stress provides the curvature in the stress-strain response. The 
cyclic hardening characteristic of the back stress is the 
consequence of an evolving limit stress. A small value of limit 
stress enables only a small amount of back stress at kinematic 
saturation. This results in a stress-strain response with sharper 
curvature, as observed in the initial, monotonic loading curve. 
The extent of separation between the applied stress and the 
back stress (i.e., the magnitude of effective stress) is a measure 
of strength associated with the drag stress. The increasing 
separation between these stresses with deformation is the 
consequence of an evolving drag stress. 

Applications 

Theoretical correlations are presented in the previous 
section; theoretical predictions are presented in this section. 
In particular, theoretical predictions are presented for 
thermomechanical and biaxial nonproportional tests for which 
experimental data exist in the literature. 

A comparison between experiment and theory is presented 
in figure 6 for an inphase thermomechanical test of copper. 
The model’s prediction is satisfactory. It does well at the tips 
of the loop, but not as well in the region of the Bauschinger 
effect. An earlier attempt to predict this experiment by the 
author (ref. 40) was less than satisfactory, because that model 
overpredicted the cold end response by some 50 percent. The 
difference is that in this model the limit stress evolves as a 
state variable, whereas in the previous model the limit stress 
varied only with temperature. Predictions were made with this 
model that both included and excluded the anisothermal T term 
in the evolutionary equations (eqs. (20) to (22)). For this 
particular application, the influence of this term on the 
predicted stress-strain response was negligible. 

The experimental stress-strain curve presented in figure 6 
has been smoothed. In reality (see ref. 40), strain avalanches 
occurred periodically along the response curve; a phenomenon 
that seems to accompany slower strain-rate experiments 
(e ^ 1x10 ~ 5 /s) done on copper. Curiously, these strain 
avalanches were nonrandom events. After cyclic saturation 
was attained, they would repeat themselves in number, degree, 
and location along the response curve from cycle to cycle. No 
attempt is made by continuum viscoplasticity to model strain 
avalanching. 

Experimental data from a cyclic, nonproportional, axial- 
torsional test done on an isotropically saturated copper 
specimen by Lamba and Sidebottom (ref. 50) are presented 
in figure 7. Theoretical predictions of this experiment are 
presented in figure 8 for three different models. One is the 
thermoviscoplastic model developed herein. The other two are 
plasticity models considered by Lamba and Sidebottom. In 
particular, they are a Krieg-Dafalias (refs. 51 and 52) two- 
surface plasticity model with a Tresca yield condition, and a 
Prager (ref. 34) plasticity model with a von Mises yield 
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AXIAL STRAIN, £ 



MECHANICAL STRAIN, £ - a AT 

(a) Experimental data (ref. 40). 

(b) Theoretical prediction. 

Figure 6.— Uniaxial, inphase, thermomechanical, stress-strain response 
for copper. Temperature, T, 200 to 500 °C; temperature rate, T, 
constant; strain rate, €, 1.5xl0' 5 + a77s. 



-200 -100 0 100 200 
AXIAL STRESS, 0, MPa 

(a) Strain control program. 

(b) Experimental stress response. 

Figure 7.— Cyclic, nonproportional, axial-torsional test of copper at room 
temperature. Test specimen cycled to isotropic saturation prior to testing. 
Path segments sequenced from 0 — 1^2 — ... ^8. Each path 
segment applied in 60 s. Data are as presented in reference 50. 





(a) Thermoviscoplastic model. 

(b) Krieg-Dafalias, two-surface, plasticity model. Data are as presented 
in reference 50. 

(c) Prager plasticity model. Data are as presented in reference 50. 

Figure 8 — Theoretical predictions of three models for the cyclic, 
nonproportional, axial-torsional test of figure 7. 


condition. 4 The viscoplastic model has a von Mises-like yield 
condition because it is a second invariant formulation. All three 
of these models were formulated such that they exhibit a strain 
hardening response that exponentially decays with continued 
directional deformation. Under proportional loading condi- 
tions, these three models produce nearly equivalent stress- 
strain responses (the difference being the rate dependence of 
the viscoplastic model versus the rate independence of the 
plasticity models). 

The viscoplastic model and the Krieg-Dafalias model both 
do a good job of representing the experimental data. The 
Prager model, however, does poorly. In addition, the Prager 
model predicts a ratcheting response that was not observed 
experimentally. The other two models do not. The Prager 
model differs from the other two models in that it does not 
incorporate the concept of limit plasticity. In particular, it does 
not require the inelastic strain rate e p to be coaxial with the 
deviatoric stress S and the back stress B at kinematic 
saturation. The Krieg-Dafalias model introduces the concept 
of limit plasticity through a second yield surface (after the 
manner of Mroz (ref. 54)) where the plastic modulus is 
zero valued; it is the limit surface. The viscoplastic model 
introduces the concept of limit plasticity through the evolu- 
tionary equation for back stress (eq. (20)) as illustrated in 
figure 1 . Here the limit stress L represents the limiting state 
of back stress. Whenever the back stress attains this value, 
the response becomes one of perfect, or limit, plasticity, 
provided that the isotropic state variables have saturated, too. 
Lamba and Sidebottom (ref. 50) conclude by stating that “For 
a given path in strain space, the path in stress space is 
determined more by limit plasticity than by any other 
consideration.” Therefore, viscopiastic models that do not 
embody the concept of limit plasticity (e.g., Bailey -Orowan 
formulations for the evolution of back stress) cannot properly 
model the plastic behavior of metals. 

Comparing the viscopiastic and Krieg-Dafalias model 
predictions with the experimental data, one notices that the 
viscopiastic model does slightly better in predicting the axial 
response, while the Krieg-Dafalias model does slightly better 
in predicting the shear response. The reason why the Krieg- 
Dafalias model predicts shear better than the viscopiastic model 
has to do with the choice of yield and limit conditions. The 
Krieg-Dafalias model of Lamba and Sidebottom used Tresca 
yield and limit conditions, while the viscopiastic model uses 
von Mises-like yield and limit conditions. A better shear 
response could be predicted for these copper data if the third 
invariant were also considered in the formulation of the 
viscopiastic model. The approach of Robinson and Ellis (ref. 
55) for incorporating the third invariant into a viscopiastic 


4 Lamba and Sidebottom (ref. 50) also considered a Ziegler (ref. 53) 
plasticity model with a Tresca yield condition. These results are not reproduced 
here because they are qualitatively equivalent to those obtained from the Prager 
model. 
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model is worthy of consideration. In contrast with the data 
of Lamba and Sidebottom (ref. 50), the copper data of Taylor 
and Quinney (ref. 56) support the von Mises yield condition 
over that of Tresca. One is therefore justified in considering 
the influence of the third invariant of stress on the inelastic 
flow of copper as a higher order effect. 

Summary of Results 

A primitive set of constitutive equations are considered from 
which a thermoviscoplastic model is developed. These 
equations assume that the material is initially isotropic, that 
the infinitesimal strain is additively decomposed into elastic 
and inelastic contributions, that the inelastic response is 
incompressible, that Onsager relations apply to the internal 
state variables, and that the inelastic strain and internal state 
variables evolve as functions of state. Associated with this 
set of constitutive equations is a thermodynamic constraint 
of intrinsic, or material, dissipativity that must always be 
satisfied. 

A thermoviscoplastic model is developed that is a special 
case of this primitive set of constitutive equations. A second 
invariant formulation with three internal state variables is 
considered. One internal variable (the back stress) accounts 
for kinematic or flow-induced anisotropic effects, while the 
other two variables account for isotropic effects. Of the two 


isotropic variables, the drag stress is a measure of yield 
strength, while the limit stress is a measure of limit strength. 
Each evolutionary equation is composed of two terms; a 
hardening term and a recovery term (either dynamic or 
thermal, but not both). The fact that there is no coupling 
between dynamic and thermal recovery terms in any evolu- 
tionary equation is a simplification that enables the model to 
reduce to both plasticity and creep models in their respective 
limits. 

There are three material functions in the model that must 
be characterized before it can be used for applications. They 
have been determined phenomenologically for copper. These 
functional forms should apply to the class of materials con- 
taining polycrystalline metals and Class II solid solution alloys 
where the effects of solute drag on strength are negligible. 

A variety of model correlations and predictions are pre- 
sented for copper under isothermal, thermomechanical, and 
nonproportional, cyclic conditions. All do a good job of 
describing the experimental data. One can therefore feel 
comfortable about using this model in real world structural 
design problems, such as regeneratively cooled combustion 
liners of rocket engines. 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, August 23, 1988 
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Appendix - Symbols 


A coefficient for inelastic strain rate 

thermodynamic forces (or affinities) 

B back stress (deviatoric) 

C power-law breakdown stress 

D drag stress 

E Young modulus 

e deviatoric strain 

F,G flow functions 

H b hardening modulus for back stress 

H b nonnegative, fourth -order, anisotropic tensor for the 
hardening modulus for back stress 
H d hardening modulus for drag stress 

H t j Onsager hardening moduli 

H ( hardening modulus for limit stress 

7 unit tensor 

k Boltsmann constant 

L limit stress 

M nonnegative, fourth-order, plastic, anisotropic tensor 

n power-law creep exponent 

Q activation energy for self-diffusion 

R thermal recovery function 

S deviatoric stress 

T absolute temperature 

T m absolute melting temperature 

Z Zener-Hollomon parameter 


a coefficient of thermal expansion 

a, thermodynamic displacements (or fluxes) 

f$ thermodynamic flux conjugate to back stress 
(deviatoric) 

5 thermodynamic flux conjugate to drag stress 
€ infinitesimal strain 

0 thermal diffusivity 

X thermodynamic flux conjugate to limit stress 
v Poisson ratio 

a Cauchy stress 

E effective stress associated with plastic flow 
(deviatoric) 

Q potential function 

Subscripts: 

i,j,k indices 

o value associated with initial or virgin states 
sat saturated condition 

a\v steady state 

2 magnitude of second-order tensor normalized for 
shear 

Superscripts: 

e elastic 

p plastic (inelastic) 

time rate of change 
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